DEFINITE INTEGRALS BY THE METHOD OF BRACKETS. 

PART 1 



IVAN GONZALEZ AND VICTOR H. MOLL 



Abstract. A new heuristic method for the evaluation of definite integrals is 
presented. This method of brackets has its origin in methods developed for 
the evaluation of Feynman diagrams. We describe the operational rules and 
illustrate the method with several examples. The method of brackets reduces 
the evaluation of a large class of definite integrals to the solution of a linear 
system of equations. 



1. Introduction 

The problem of analytic evaluations of definite integrals has been of interest to 
scientists since Integral Calculus was developed. The central question can be stated 
vaguely as follows: 

given a class of functions 5 and an interval [a, 6] C M, express the integral of f ^ ^ 

I = f{x)dx, 



in terms of the special values of functions in an enlarged class ©. 

For instance, by elementary arguments it is possible to show that if is the 
class of rational functions, then the enlarged class © can be obtained by including 
logarithms and inverse trigonometric functions. G. Cherry has discussed in [17], 
[18] and [19] extensions of this classical paradigm. The following results illustrate 
the idea: 

/i;^4«(-*-^^" + "4'"--'>' 

but 
(1.2) 



ln(a;2 _ 1) 

can not be written in terms of elementary functions and the logarithmic integral 

X dx 
Inx 

that appears in (1.1). The reader will find in [16] the complete theory behind 
integration in terms of elementary functions. 



(1.3) liix) 
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Methods for the evaluation of definite integrals were also developed since the 
early stages of Integral Calculus. Unfortunately, these are mostly ad-hoc procedures 
and a general theory needs to be developed. The method proposed in this paper 
represents a new addition to these procedures. 

The evaluations of definite integrals have been collected in tables of integrals. 
The earliest volume available to the authors is [8], compiled by Bierens de Haan 
who also presented in [9] a survey of the methods employed in the verification of 
the entries. These tables form the main source for the popular volume by I. S. 
Gradshteyn and I. M. Ryzhik [39]. Naturally any document containing a large 
number of entries, such as the table [39] or the encyclopedic treatise [51], is likely 
to contain errors. For instance, the appealing integral 

(1.4) ^ ' 



POO 

Jo n. 



1/2 2^/6' 



(1 + a;2)3/2 ^^(x) + 
with 

4a;^ 

(1-5) ^(-) = i + 3(TT^' 

that appears as entry 3.248.5 in [38], the sixth edition of the mentioned table, is 
incorrect. The numerical value of I is 0.666377 and the right hand side of (1.4) is 
about 0.641275. The table [39] is in the process of being revised. After we informed 
the editors of the error in 3.248.5, it was taken out. There is no entry 3.248.5 in 
[39] . At the present time, we are unable to evaluate the integral I. 

The revision of integral tables is nothing new. C. F. Lindman [47] compiled 
a long hst of errors from the table by Bierens de Haan [10]. The editors of [39] 
maintain the webpage 

http : //www . mathtable . com/gr/ 

where the corrections to the table arc stored. The second author has began in 
[2, 3, 55, 56, 57, 58, 59, 60] a systematic verification of the entries in [39]. It is in 
this task that the method proposed in the present article becomes a valuable tool. 

The method of brackets presented here, even though it is heuristic and still lacking 
a rigorous description, is quite powerful. Moreover, it is quite simple to work with: 
the evaluation of a definite integral is reduced to solving a linear system of equations. 
Many of the (nitries of [39] can be derived using this method. The basic idea behind 
it is the assignement of a bracket (a) to any parameter a. This is a symbol associated 
to the divergent integral 



/.oo 

(1.6) / x'^-^dx. 

Jo 



The formal rules for operating with these brackets are described in Section 3 and 

their justification is work-in-progress, we expect to report in the near future. The 
rest of the paper provides a list of examples illustrating the new technique. 
Given a formal sum 

oo 

(1.7) /(x) = ^a„x«"+'5-i 

n=0 
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we associate to the integral of / a bracket series written as 

(1.8) / f{x)dx = y2an{an + f3), 

Jo 

to keep in mind the formahty of the method described in this paper. Convergence 
issues are ignored at the present time. Moreover only integrals over the half-line 
[0,00) will be considered. 

Note. In the evaluation of these formal sums, the index n € N will be replaced by 
a number n* defined by the vanishing of the bracket. Observe that it is possible 

that n* e C. For book-keeping purposes, specially in cases with many indices, we 

00 

write instead of the usual After the brackets are eliminated, those indices 

n n=0 

that remain recover their original nature. 

The rules of operation described below assigns a value to the bracket series. The 
claim is that for a large class of integrands, including all the examples described 
here, this formal procedure provides the actual value of the integral. Many of the 
examples involve the hypergeometric function 

F ._ ("0" (°'2)n • • • (ap)n ^" 

This series converges absolutely for all 2 e C if p < (7 and for \z\ < 1 \i p = q + 1. 

The series diverges for all z 7^ if p > g + 1 unless the scries terminates. The 
special case p = g + 1 is of great interest. In this special case and with \z\ = 1, the 
series 

(1.10) g+iFq(oi,... ,ag+i;6i,... ,hq\z) 

converges absolutely if Re (^ 6^ — ^ Oj) > 0. The series converges condition- 
ally if z = e** 7^ 1 and > Re(^6j — X^Oj) > —1 and the series diverges if 
Re(E^.-Eai)<-l- 




Figure 1 . The triangle 



The last section of this paper employs the method of brackets to evaluate certain 
definite integrals associated to a Feynman diagram. From the present point of view, 
a Feynman diagram is simply a generic graph G that contains E + 1 external lines 
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and N internal lines or propagators and L loops. All but one of these external lines 
are assumed independent. The internal and external lines represent particles that 
transfer momentum among the vertices of the diagram. Each of these particles 
carries a mass rrii > for i = 1, N . The vertices represent the interaction of 
these particles and conservation of momentum at each vertex assigns the momentum 
corresponding to the internal lines. A Feynman diagram has an associated integral 
given by the parametrization of the diagram. For example, in Figure 1 we have 
three external lines represented by the momentum Pi, P2, P3 and one loop. The 
parameters aj are arbitrary real numbers. The integral associated to this diagram 
is given by 

r(ai)r(a2)r(a3) X Jo Jo {xi+X2+X3)''/^ 

( 2 , 2 , 2\ ( CiiPl + 2C12P1 ■ Pi + C22Pi\ , 

expfxiTO-, + X2mr, + xsmo) exp dx, 

V XI+X2+XS J 

where dx := dxi dx2 dx^. The evaluation of this integral in terms of the variables 
Pi & M.'^ , rui & M. and Oj G M is the solution of the Feynman diagram. The functions 
Cij are polynomials described in Section 13. 

The method of brackets presented here has its origin in quantum field theory 
(QFT). A version of the method of brackets was developed to address one of the 
fundamental questions in QFT: the evaluation of loop integrals arising from Feyn- 
man diagrams. As described above, these arc directed graphs depicting the interac- 
tion of particles in the model. The loop integrals depend on the dimension D and 
one of the (many) intrinsic difficulties is related to their divergence at £> = 4, the 
dimension of the physical world. A correction to this problem is obtained by taking 
D = A — 2e and considering a Laurent expansion in powers of e. This is called the 
dimensional regularization [11] and the parameter e is the dimensional regulator. 

The method of brackets discussed in this paper is based on previous results by 
I. G. Halliday, R. M. Ricotta and G. V. Dunne [26], [27] and [40]. The work in- 
volves an analytic extension of D to negative values, so the method was labelled 
NDIM (negative dimensional integration method). The validity of this continua- 
tion is based on the observation that the objects associated to a Feynman diagram 
(loop integrals as well as the functions linked to propagators) are analytic in the 
dimension D. A. Suzuki and A. Schmidt employed this technique to the evaluation 
of diagrams with two loops [68], [69]; three loops [71]; tensorial integrals [70] and 
massive integrals with one loop [66], [72], [67]. An extensive use of this method as 
well as an analysis of the solutions was provided by C. Anastasiou and E. Glover 
in [5] and [6]. The conclusion of these studies is that the NDIM method is inade- 
quate to the evaluation of Feynman diagrams with an arbitrary number of loops. 
The proposed solutions involve hypergeometric functions with a large number of 
parameters. By establishing new procedural rules I. Gonzalez and I.Schmidt [36] 
and [37] have concluded that the modification of the previous procedures permits 
now the evaluation of more complex Feynman diagrams. One of the results of [36], 
[37] is the justification of the method of brackets in terms of arguments derived 
from fractional calculus. The authors have given NDIM the alternative name IBFE 
(Integration by Fractional Expansion). 

From the mathematical point of view, the NDIM method has been used to pro- 
vide evaluation of a very limited type of integrals [64], [65]. The examples presented 
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in this paper show great flexibility of the method of brackets. A systematic study 
of integrals arising from Feynman diagrams is in preparation. 

2. A DETOUR ON DEFINITE INTEGRALS 

The literature contains a large variety of techniques for the evaluation of definite 
integrals. Elementary techniques arc surveyed in classical texts such as [28] and 
[32] . The text [7] contains an excellent collection of problems solved by the method 
of contour integration. The reader will find in [13] a discussion of several elementary 
analytic methods involved in the evaluation of integrals. 

It is hard to predict the type of techniques required for the evaluation of a specific 
definite integral. For instance, [73] contains a detailed account of the proof of 

(2.1) / lnlntana;da;=-lnN^ — TTT^ ' 

that appears as formula 4.229.7 in [39]. This particular example involves the use of 
X-functions 

(2-2) L^is) ^±^, 

n=l 

where X is a character. This is a generalization of the Riemann zeta function ^(s) 
(corresponding to x = !)• Vardi's technique has been extended in [52] to provide a 
systematic study of integrals of the form 

(2.3) Iq= [ Q{x)lnlnl/xdx, 

Jo 

that gives evaluations such as 

(2.4) I (. + ;). = 2(-l-'2 + ,-ln. + ln2), 

f°° TTT^ 3 7r^ln2 

(2.5) / Inxlntanhxrfa; = tC'(2) H 777—. 

Jo 8 4 12 

Here 7 = — r'(l) is Euler's constant. 

A second class of examples appeared during the evaluation of definite integrals 
related to the Hurwitz zeta function 

00 

(2-6) C{z,q)=y- -. 

In [29] the authors found an evaluation that generalizes the classical integral 



/ lnT{q)dq = \nV2^, 
Jo 



(2.7) 
namely 

where Li is in (2.7) and A = 7+ln 27r. The natural next step, namely the evaluation 
of 



(2.9) L3:= / ln"r(g)rfg. 



Jo 
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remains to be completed. In [30] and [31] the reader will find a relation between 
Ls and the Tornheim sums T(m, k, n), for m,k,n G N. These sums are defined by 



oo oo 



(2.10) T(a,b,c) := > / . , 
The special case 

(2.11) T{n,0,m)= J2 



corresponds to the multiple zeta value (MZV) ({n,m) of depth 2. The MZV is 
given by 

(2-12) ({ni,n2,--- ,nr) := V :^7Jr77J7^ 

h>h>->jr^^ -^2 

where the parameter r is called the depth of the sum. These series were initially 
considered by Euler and have recently appeared in many different places. The 
reader will find in [44] a description of how these sums are connected to knots and 
Feynman diagrams. These diagrams arc a very rich source of interesting integrals. 
The last section of this paper is dedicated to the evaluation of some of these integrals 
by the method of brackets. 

The computation of hyperbolic volumes of 3-manifolds provides a different source 
of interesting integrals. Mostow's rigidity theorem states that a finite volume 3- 
manifold has a unique hyperbolic structure. In particular its volume is a topological 
invariant. An interesting class of such 3-manifolds is provided by hyperbolic knots 
or link complement in S^. The reader will find information about this topic in the 
articles by C. Adams and J. Weeks in [53]. It turns out that their hyperbolic struc- 
ture can be given in terms of hyperbolic tetrahedra [1] . Milnor [54] describes how 
the volume of these tetrahedra can be expressed in terms of the Clausen function 

u 

(2.13) Ch{0):=- log]2sin-|du. 

Jo 2 

The reader will find in [49] a discussion on arithmetic properties of 3-manifolds. In 
particular, Chapter 11 has up to date information on their volumes. 

Zagier [74] provided an arithmetic version of these computations in his study of 
the Dedekind zeta function 

for a number field K that is not totally real. Here N{a) is the norm of the ideal a 
and the sum runs over all the nonzero integral ideals of K. In the case of totally 
real number fields a classical result of Siegel shows that Cx(2m) is a rational mul- 
tiple of 7r2"™/-\/!D, where n and D denote the degree and the discriminant of K, 
respectively. Little is known in the non-totally real situation. Zagier [74] proves 
that Cif (2) is given by a finite sum of values of 

(2.15) A(.):=£^log^.. 
The function A can be written as 

(2.16) A{x) = Cl2(7r-2tan-^a;)-Cl2(7r). 
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Morever, he conjectured that Cfe(2TO) can be given in terms of 

-)2m—l r°° j.2m—l ^j. 



22m— 1 ro 



X sinh t + x^'^ cosh t 



The conjecture is estabUshed in the special case where K is an abelian extension of 
Q. The example K = Q{^/^) yields 



7r2 



(2.18) Cq(V=7)(2) = ^ (cot ^)+A (cot '-f)+A (cot f )) 
and also 

(2.19) Cq(V^) (2) = ^ (2A{V7) + AiV7 + 2^3) + ^(^7 - 2x/3)) , 
leading to the new Claussen identity 

A (cot f )+A (cot ^)+A (cot ^) = f (2A{V7) + A{V7 + 2\/3) + A{V7 - 2\/3)) . 

Zagier stated in 1986 that there was no direct proof of this identity. To this day this 
has elluded considerable effort. The famous text of Lewin [46] has such parametric 
identities but it misses this one. R. Crandall [23] has worked out a theory in which 
certain Claussen identities are seen to be equivalent to the vanishing of log-rational 
integrals. 

J. Borwein and D. Broadhurst [15] identified a large number of finite volume 
hyperbolic 3-manifolds whose volumes are expressed in the form 

Here i^T is a field associated to the manifold 371 (the so-called invariant trace field) 
and n and D arc the degree and discriminant of K, respectively. The authors offer 
a systematic numerical study of the rational numbers — . The identity of Zagier 
described above yields the remarkable identity 



r^'ln ^5^ii±^ dt = A{V7) + Ui^+'^^) + Ui^-'^^)- 
Jtt/3 tanf-v7 2 2 



tan t + \/l 



This example corresponds to the link 6f with discriminant D = —1 . Zagier's result 
gives ^ = 2 in (2.20). 

Coffey [20], [21] has studied the integral above, that also appears in the reduction 
of a multidimensional Feynman integral [48] . The goal is to produce a more direct 
proof of Zagier remarkable identity as well as the many others that have been 
numerically verified in [15]. 

The subject of evaluation of definite integrals has a rich history. We expect that 
the method of brackets developed in this paper will expand the class of integrals 
that can be expressed in analytic form. 
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3. The method of brackets 

The method of brackets discussed in this paper is based on the assignment of 
a bracket (a) the parameter a. In the examples presented here a € R, but the 
extension to a € C is direct. The formal rules for operating with these brackets are 
described next. 

Definition 3.1. Let / be a formal power series 

oo 

(3.1) f{x) = J2anX"^^^-K 

n=0 

The symbol 

/•oo 

(3.2) / f{x)dx = y2an{an + p) 

represents a bracket series assignement to the integral on the left. Rule 3.2 describes 
how to evaluate this series. 

Definition 3.2. The symbol 

(-1)" 

(3.3) 0„ := ^^^^ 
will be called the indicator ofn. 

The symbol (pn gives a simpler form for the bracket series associated to an integral. 
For example, 



(3.4) 



/■oo 

/ x°'~^e~'' dx = 'S2(t>n{n + a). 
Jo „ 



The integral is the gamma function T(a) and the right-hand side its bracket expan- 
sion. 

Rule 3.1. For a G C, the expression 

(3.5) {ai+a2 + --- +ar)" 

is assigned the bracket series 



(-a + mi H h rrir) 



(3.6) J2 ^i,2,-,rar---a^ 

mi,--- ,mr- ^ 

where (j)i,2,---,r is a short-hand notation for the product (?i'mi< 

Rule 3.2. The series of brackets 

(3.7) ^(/)„/(n)(on + 6) 

n 

is given the value 

(3.8) i/(n*)r(-n*) 
where n* solves the equation an + 6 = 0. 



'1712 H^m-r • 
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Rule 3.3. A two-dimensional series of brackets 

(3.9) ^ 0ni,n2/("l:"2)(aiini + 012^2 + Ci)(a2ini +022^2 + C2) 

is assigned the value 

(3.10) \ finl,n*,)Ti-nl)T{-n*^) 

where (n*,n2) is the unique solution to the linear system 

(3.11) aiiTii + ai2n2 + ci = 0, 

a2ini + a22n2 + C2 = 0, 

obtained by the vanishing of the expressions in the brackets. A similar rule applies 
to higher dimensional series, that is, 

• • • ,r/(^l) • ■ ■ ) ^r)(ttlini + • • • airUr + Ci) • • • {ttririi + ■ ■ ■ arrlT-r + C^) 

is assigned the value 

(3.12) ^i^-^^'^i' ■ ■ ■ ' OI'^-O • • • 

where A is the matrix of coefficients (a^) and {n* } is the solution of the linear 
system obtained by the vanishing of the brackets. The value is not defined if the 

matrix A is not invertible. 

Rule 3.4. In the case where the assignment leaves free parameters, any divergent 
series in these parameters is discarded. In case several choices of free parameters 
are available, the series that converge in a common region are added to contribute 
to the integral. 

A typical place to apply Rule 3.4 is where the hypergeometric functions pFg, with 
p = q + 1, appear. In this case the convergence of the series imposes restrictions 
on the internal parameters of the problem. Example 13.2, dealing with a Feynman 
diagram with a bubble, illustrates the latter part of this rule. 

Note. To motivate Rule 3.1 start with the identity 

1 1 /.oo 

(3.13) _ = -_/ x"-'^e-^'' dx, 

and apply it to A = ai H + a^. to produce 

1 f°° 

(ai H 1- Or)" = TT, T / a;~'*~-^exp [-(fli H \- ar)x] dx 

r(-a) Jo 



r( 

Expanding the exponentials we obtain 



1 

/ x-'^-^e-"'"' ■■■e-"-'' dx. 

-a) Jo 



1 /-oo 

(ai + • • • + = T y---y ^h- ,rar •••<'■/ a;-"+"i+-+™'-i dx 

and thus 

(3.14) (a, + ... + a.)«^^...^c^ (-a + + • • • + m.) 



T{-a) 
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This is Rule 3.1. 



4. Wallis' formula 

The evaluation 

J _ r _ ^ ^2m^ 

^ ' Jo (1 + ar2)'"+i - 22^+1 V m , 

is historically one of the earliest closed-form expressions for a definite integral. The 
change of variables x = tan 9 converts it into its trigonometric form 

(4.2) J,,^:=£'\os^-9de=^^^"' 
An elementary argument shows that J2,m satisfies the recurrence 

(4.3) J2 TO = n •^2,TO-l 

2m 

and then one simply checks that the right hand side of (4.2) satisfies the same 
recurrence with matching initial conditions. A second elementary proof of (4.1) is 
presented in [14]: using cos^ 9 = |(1 + cos 20) one obtains the recurrence 



Lm/2J 

m 



(4-4) J2,TO=2- y: qj^ 

and the inductive proof follows from the identity 

Lm/2J 

(4.5) J2 2 



2,i, 



i=Q 



Jm\ /2A ^ 2-™P"^ 
2i)\i J \m J' 



This can be established using automatic methods developed by H. Wilf and D. 
Zeilberger in [61]. 

The proof of Wallis' formula by the method of brackets starts with the expansion 
of the integrand as 



(4.6) (l + a;2)-™-i4^^<^i,2 



m + l + ni+n2) ^2n2^ 



r(m+l) 

The corresponding integral J2,to is assigned the bracket series 

(4.7) J2,TO =VV</>i,2 7^7-^— T(TO+l + ni +n2) (2712 + 1). 

^ — ' ^ — ' 1 (m + 1) 

Rule 3.2 then shows that 

T - i r(-nt)r(~n*) 
^^■^> ■^''™-2 r(m+l) ' 

where (n^, 712) is the solution to the linear system of equations 

(4.9) m + 1 + ni + 712 = 0, 

2712 + 1 = 0. 
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Therefore nl = — (to + ^) and rij = We conclude that 

(4.10) J2,m = 



r(TO + i)r(i) 



2r(m) 

This is exactly the right-hand side of (4.1). 

5. The integral representation of the gamma function 
The exponential in the integral 

(5.1) / = / x'^-^e-'^dx 

Jo 

is expanded in power series to obtain 

(5.2) x-'e-^ = J2 ^-^^^^ = E 

Therefore, the integral (5.1) gets assigned the bracket series 

(5.3) 7 = 5^0„(a + n). 

n 

Rule 3.2 assigns the value T{a) to (5.3). This is precisely the value of the integral: 

(5.4) / x^'-^e-'' dx = T{a). 

Jo 

Rule 3.2 was developed from this example. 

6. A Fresnel integral 
In this section we verify the evaluation of Fresnel integral 

(6.1) / sin(aa; ) dx 



2\/2a 

The reader will find in [7] the standard evaluation using contour integrals and other 
elementary proofs in [33] and [45]. 

In order to apply the method of brackets, use the hypergeometric representation 



sin^ _ 

— O-C^l 



that can be written as 



(6.2) sin z = '^(j)n 



■ 3. _£l 

' 2' 4 



n=0 \2Jn 

Therefore 

j2n+l 

(6.3) / siii{ax^) dx = y (l)njj^ (4n + 3). 



/ sin(aa;2 ) da; = ^ /gx ■„ 

■^0 „ \2)n^ 



According to Rule 3.2, the assignment of the right-hand side is obtained by evalu- 
ating the function 
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at the solution of 4n* +3 = 0. Therefore the integral (6.1) has the value 

where the factor j comes from the term 4n + 3 in the bracket. Using {a)m = 
r(a + m)/r(a), we obtain 



(6.5) 



3\ 2r(|) 



-3/4 



Wo conclude that the assigned value is 7r/2v2a. As expected, this is consistent 
with (6.1). 

The method also give the evaluation of 

/■OO 

(6.6) / = / x''-^sm{ax'=)dx. 

Jo 

The change of variables t = transforms (6.6) into 

1 

(6.7) I=- t''/''-^sm{at)dt, 

c Jo 

and this is formula 3.761.4 in [39] with value 

(6.8) ^°°.^-sin(a.^)c^.= M,ing). 

To verify this result by the method of brackets, start with the expansion 

CO 2n+l 

(6.9) x"-' Max^) = ^ 

71=0 \2)n 

and associate to it the bracket series 

(6.10) / x^~'^sm{ax'')dx = ^4>n-7jT — ^(2nc + c + 6). 

•^0 n (sin 2" 

Apply Rule 3.2 to obtain 

1 a^n.+i 

y ' 2c (i) . 22"* ^ " 

where n* solve 2nc + 6 + c = 0; that is, n* = -1/2 - 6/2c. Then (6.11) yields 

(Q12) I- ^^^^^''^'^ r(i + A) 

' caVcr(l- A)' ^2 + 2cj- 

with X = b/2c to transform (6.12) into (6.8). To transform (6.12) into (6.8), simplify 

(6.12) using the reflection formula 

(6.13) r(x)T(l-x) = 

sin wx 

and the duplication formula 

(«») = 

with X = b/2c. 



EVALUATION OF INTEGRALS BY BRACKETS 



13 



Note. The method developed by Flanders [33] is based on showing that 

POO POC 

(6.15) F{t):= e-'^'cosa^^da; and G(i) := / e"*^' sinx^ 

Jo Jo 



10 

satisfy the functional equation 
(6.16) F^{t) - G^{t) = 2F{t)G{t) = 



4(1 +i2)- 

The latter can be solved to obtain the values 



(6.17) F{t) = J-\ r ^ , . andG(t) 



1 + i2 ' V 8 V 1 + ' 

A second elementary proof was obtained by Leonard [45]. Converting (6.15) into 
the Laplace transform of cosa;/2y^ and sinx/2y^ respectively, he shows that 

The evaluation of these integrals described in [45] , is elementary but long. A shorter 
argument follows from the formula 

(6.19) ' 



a) := / 
Jo 



c4 + 2ax2 + l 2^2(1 + a)' 
Indeed, the values 



G{t) = (1 + e)-^'\f y-^^^) and F{t) = (1 + t'r^l^Git/ ^/l^) + tG(t), 

follow from (6.19) by a change of variable v = {1 + t^Y^'^u. The evaluation of the 
quartic integral (6.19) by the method of brackets is discussed in detail in Section 
12. 

7. An INTEGRAL OF BETA TYPE 

In this section we present the evaluation of 

r°° dx 
~ Jo {E + Fx^y 

The change of variables x = C^l^t^l^ , with C = E/F, yields 
(7.2) ^ C [ t dt 



(7.1) 



hE- Jo (1 + 

where u= {a + l)/b. The new integral evaluates as B{c — u, u) where B{x, y) is the 
classical beta function; see [39], formula 8.380.3. We conclude that 

(7.3) i=^^B{c-u,u). 

To evaluate this integral by the method of brackets, the integrand {E + Fx'')~'' 
is expanded as 

ni pn2 3»2 (c + ni +n2) 



(7.4) ^^<^i,2i;"^F"=a;' 



r(c) 
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Replacing in (7.1) wc obtain 

I L yy^i,E"^F"- (c + ni+na) /•=^^,+,„,+i_i 
ft' ' ^^"^ 
= E E '^i>2^"'^"' + "1 + "2)(a + bn2 + 1). 

To obtain the value assigned to the two dimensional sum, solve 

c + ni + n2 = 
a + 6n2 + 1 = 0, 
to produce the solution n\ = — c and rv^ = — . Therefore 

(7.5) / = ^^E<F^'^T{-nl)T{-n*^), 

and this reduces to the value in (7.3). 

8. A COMBINATION OF POWERS AND EXPONENTIALS 

In this section we employ the method of brackets and evaluate the integral 

x°'~^ dx 



-I 

Jo 



with a, (3, 7, A, B, C G R. To evaluate this integral we consider the bracket series 
(8.2) {A + BeMCx^))'^ = ^"^■B"^exp(Cn2X^) ^ | ""'^ . 

ni ,n2 

The exponential function is expanded as 

713=0 ^ " ' 



X:C"3(-n2) 



Therefore, the integral (8.1) is assigned the bracket scries 

^ni C-ns (_n2)"3 (a + /Jns) (7 + ni + 712) 



^1,2,3" 



r(7) 

ni.n2,"3 

The vanishing of the two brackets leads to the system 

a + /3n3 = 
7 + ni + n2 = 0, 

and we have to choose a free parameter between m and n2- Observe that na = 
—a/ (3 is determined by the method. 



Choice 1: take n2 to be free. Then nl = —7 — 712 and = —a/ (3. This leads to 

_S"=r(a//3)r(7 + n2) 
< A' 

n2- 



(8.3) I=J2 



A7+n2Ca//3/3r(7)(-n2)«/^ ' 
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This is impossible due to the presence of the term '''' leading to a divergent series. 
These divergent series are discarded. 

Choice 2: take rii as the free variable. Then = —a/P and = —7 — ni. This 
time we obtain 



(8.4) I = 



r(7) B7C«//3/3 ' r(l + ni) (7 + ni)"//3' 

m— ^ ' 



This formula cannot be expressed in term of more elementary special functions. 
In the special case 7 = 1 we obtain 

(8.5) / = -IMpolyLog(i/, -a/b). 

with u = a/ (3. The polylogarithm function appearing here is defined by 

(8.6) PolyLog(2,fc):=5^— . 

n— 1 

Specializing to A = B = C = a = ^ — \ and /? — 2 we obtain 

Of course, this integral can be evaluated by simply expanding the integrand as a 
geometric series. 

9. The Mellin transform of a quadratic exponential 
The Mellin transform of a function /(x) is defined by 

/■CX) 

(9.1) M{f){s):= x'-'f{x)dx. 

Jo 

Many of the integrals appearing in [39] are of this type. For example, 3.462.1 states 

that 

(9.2) 

M (e-'^^-'-^-) (s) = x'-'e-('-"-"'-dx = (2/3)-^/2r(s)e^'/(««i?_, (^-^^ . 
Here Dp{z) is the parabolic cylinder function defined by (formula 9.240 in [39]) 



r((i-p)/2)^ ^ V 2' 2' 2; r(-p/2) 

A direct application of the method of brackets gives 
(9.3) / a;"-ie-^^'-''^da; = ^^(/)i,2/3"'7"'(s + 2ni+n2). 

The equation s + 2ni + n2 = gives two choices for a free index. Taking 
= — 2ni — s leads to the series 

00 -. ^ 1 //I /3\ ^1 



T. rin!+i){ 7^0 ~ „?,r(niVi)( 7^) 



ni=0 ^ / \ / / ni=0 

's S + 1 



s + 1 



4/3 
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This choice of a free index is excluded because the resulting series diverges. The 
second choice is nj = —nijl — s/2 and this yields the series 

^^■^^ ^ r(n2 + i)^ ( Y + 2) ' 

where p = — To write (9.4) in hypergeometric form we separate it into two 
sums according to the parity of 722 and obtain 



The identity 

gives the final representation of the sum as 



This is (9.2). 



r(|)...(|,i;i.^)..r(l±£)... 



1 + s 3 1 . 



2 ' 2' 2 



The special case s = 1 gives 

(9.7) e-^-^-^- = 273 1^1 ihh hp') + p i^^i (1; I; ip')] ■ 

The first hypergeometric sum evaluates to e^^/'^^ and using the representation of 
the error function 

(9.8) erf(a;) = -= e"*' dt 
as 

(9.9) erf(a;) = ^e-'^\Fi(l;f;.T2), 

(given as 8.253.1 in [39]) we find the value of the second hypergeometric sum. The 
conclusion is that 

This can be checked directly by completing the square in the integrand. 

10. A MULTIDIMENSIONAL INTEGRAL FROM GRADSHTEYN AND RYZHIK 

The method of brackets can also be used to evaluate some multidimensional 
integrals. Consider the following integral 

,00 ,oc. ,00 „i-l^p,-l . .^p„-l ^^^^^2 ■■■dXn 



^^^■^^ 7o 7o Jo {l + {riXiY^+--- + {rnXnY-r 

which appears as 4.638.3 in [39] with an incorrect evaluation. 
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The first step in the evaluation of /„ is to expand the denominator of the inte- 
grand using Rule 3.1 as 



1 



Next the integral is assigned the value 

(10.2) 7„4 Y: 'Po,- .n f[{r,x,r'^ + ^° + . • + ^"^ nfa- + 
The evaluation of this bracket sum involves the values 

n 

(10.3) ko = -S + J2— and kj = for l<j<n. 

j=i 1j 

We conclude that 



The table [39] has the exponents of Vj written as pjqj instead oi pj. This has now 
been corrected. 

11. An EXAMPLE INVOLVING BESSEL FUNCTIONS 

The Bessel function Ju{x) is defined by the series 



oo 

(11.1) j^(:,) = ^(_i)fc. 



z 



2k+v 



2^ L^"- > 22fefc!r(!/ + fc + l)' 



and it admits the hypergeometric representation 



The method of brackets will now be employed to evaluate the integral 

/•oo 

(11.3) \ x-^J^{ax)J^,{Px)dx. 

Jo 

Three integrals of this type form Section 6.574 of [39] . 

Replacing the hypergeometric form in the integral, we have 

(f )-(!)" 



I = 



r(z/ + i)r(M + 1) 



' 4"i+»2 (i. + 1)„^ (^ + 1) 
Therefore, the bracket series associated to the integral (11.3) becomes 



I = 



EE /':' r (2ni + 2n2-A + ^ + M+l). 
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The vanishing of the brackets yields the value nl = ^(A — i^ — /x— 1) — n2 and 
it follows that 

" r(z. + l)r(M+l)„^o4»J+»^ + + 2 

Writing the Pochhammer symbol {v + in terms of the gamma function we 
obtain 



2^r(/x + +l) 

In order to write this in hypergeometric terms, we start with 



2^r(M + +i) 

. ^y,, (^(z. + /x-A+l))„,(/3Va2)"- 
,t^n (i(A + i.-M + !))-„, (/x + l)„.r(n2 + l)' 



X 

n2=0 "^2 



and use the identity 

(-1)" 

(11.4) (c)_„= ^ ' 



(1 - C)r 

to obtain 



7 = 



r(i(z/ + /x- A + 1) 



2^ r(M + i)r(i(A + !/-/i + i)) 

E (^(1 - ^ - - + /^))- (M- + - A + l)U (^+i)^^r(n2 + l) j 



that can be written as 

/J'^a^-'^-i T{^{iy + n - X + 1)) 



2^ r(^ + i)r(i(A + i^-A* + i)) 



X 2-Fl 



2 ^ 

i(l-A-i/ + At) i(z/ + /i-A + l) 

/U+ 1 



This solution is valid for |/3^/a^| < 1 and it corresponds to formula 6.574.3 in 
[39]. The tabic contains an error in this formula, the power of P is written as u 
instead of /x. To obtain a formula valid for |/3^/a^| > 1 we could proceed as before 
and obtain 6.574.1 in [39]. Alternatively exchange {i^,a) by and use the 

formula developed above. 

12. A NEW EVALUATION OF A QUARTIC INTEGRAL 

The integral 
(12.1) A^o,4(a;m) := / 

^0 



(a;4 + 2aa;2 + 1)™+! 
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is given by 

(12.2) 7Vo,4(a,m) = ; ^^^"^ 

^ [2(a+l)]'"+2 

where is the polynomial 

m 

(12.3) P„(a) =^di,„a', 
with coefficients 

«-.--|;^'(^::r)C"r)(: 

The sequence {d;,™ : < Z < m} have remarkable arithmetical and combinatorial 
properties [50]. 

The reader will find in [4] a survey of the many difTcrcnt proofs of (12.2) available 
in the literature. One of these proofs follows from the hypergeometric representation 
(12.5) 



1 1 



A^o,4(a, m) = 2"- 2 (a + \)-^-^B {2m + |) 2P1 



— m TO + 1 



TO + I 



1 - a 



New proofs of this evaluation keep on appearing. For instance, the survey [4] 
does not include the recent automatic proof by C. Koutschan and V. Levandovskyy 
[43]. The goal of this section is to provide yet another proof of the identity (12.2) 
using the method of brackets. 

The bracket series for I = No^4{a,m) is formed by the usual procedure. The 
result is 

(12.6) 1= ^ Yl <^i,2,3(2a)"=(4ni+2n2 + l)(TO + l + ni+n2+n3). 

ni,n2,n3 

The expression (12.6) contains two brackets and three indices. Therefore the 
final result will be a single series on the free index. We employ the following no- 
tation: / is the original bracket series, the symbol Ij denotes the series / after 
eliminating the index rij. Similarly lij denotes the series / after first eliminating 
rii (to produce li) and then eliminating nj. 

Case 1: 7X3 is the free index. Eliminate first rii from the bracket (4ni + 2n2 + 1) 
to obtain = —^n2 — j- The resulting bracket series is 

(12.7) /, 4 ^ ^^.3^^^5g^(m + I + In^ + n,). 

n2,n3 ^ ' 

The next step is to eliminate ni to get n% = —2m — | — 2n3 and obtain 

1 °^ 6 

^''■'^ '^'^ = 2r(TO + l)(2a)^-^3/2 E J^n-m - I - n^Wm + | + 2n3). 

In order to simplify these expressions, we employ 
(12.9) T{x + m) = {x)mT{x), T{x - m) = (-l)"T(a;)/(l - x)^ 
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and 

(12.10) {X)2m = 2 

for a; e M and m e N. We obtain 



r(-m - h - ns) 



(_l)n3r(-i -m) 
(i+TO)n3 



and 



These yield 



r(2r 



2n3) = r(2m + |)(m + |)„3(m + |)„3 22"3. 



...... . _ r(-i - m)r(2m + I) ^ (m + 3/4)„3 (m + 5/4)„3 

^ • ^ 2r(m + l)(2a)2™+3/2 (m + 3/2)„3n3! ' 



or 



r(-i -TO)r(2m + 1) 



m + 



m + 



m + 



1 



(12.12) 7i,2 2r(m + l)(2a)2™+3/2 

Note. The reader can check that /i,2 = h,!, so the value of the sum for the quartic 

integral docs not depend on the order in which the indices ni and n2 arc eliminated. 
The reader can also verify that this occurs in the next two cases described below; 
that is, /i,3 = 73,1 and 72,3 = -^3,2. 

Case 2: m is the free index. A similar argument yields 



(12.13) 



r(m + i)r(i) 

2r(m+l)(2a)V2 i 



I 1 



m 



Case 3: n2 is the free index. Eliminate rii from the bracket series (12.6) to produce 



(12.14) 



52,3- 



(2 a)"^r(|n2 + \) 
4r(m + 1) 



(m + I + in2 + n3), 



and now eliminate to obtain = —m — j — ^n2. This yields 

'''' = irwT) ^ ^^(^"^ + + 1 + 

n2=0 

In order to obtain a hypergeometric representations of these expressions, we 
separate the last series according to the parity of n2- 



h,3 



^ OO 

i V 

■71 -i- ^\ ^ 



4r(„. + 1) IS t-'"^ + + + « 

^ ' n2=0 

1 [2af 
4r(m+l)^^^(2n2 + l)! 



^2r^2+l 



— yr(n2 + |)r(n2+m+|). 



Using the standard formulas (12.9) and (12.10), we can write this in the form 



^1,3 



2-^1 



r(i)r(m + | ) ^ (\ m + 

4r(m + 1) 

ar(|)r(m + |) 

2r(m + 1) 



2-^*1 I * 3 
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In summary: we have obtained three series related to the integral Na^4{a,m). 
The series /i,2 and 12,3 are given in terms of the hypergcomctric function 2F1 with 
last argument These series converge when > 1. The remaining case /i,3 

gives 2F1 with argument a^, that is convergent when < 1. Rule 3.4 states that 
we must add the series /i_2 and 72,3 to get a valid representation for > 1. In 
conclusion, the method of brackets shows that 



No,4{(i, m) = 



r(i)r(m + |) 

4r(m + 1) 
ar(|)r(m+|) 
2T{m + 1) 

r(i)r(m + i) 



2FI 

2F1 
2FI 



m + 



+ 



m - 



+ 



2^/2aT{m + 1) 

r(-i)r(2m + |) 

2(2a)2"+3/2r(m + l) 



for < 1, 



2-^1 



m + 



4 1 for > 1. 



The continuity of these expressions at a = 1 requires the evaluation of 2-F1 (a, b;c;l) 
Recall that this is finite only when c > a + b. In our case, we have four hypcrgc- 
ometric terms and in each one of them, the corresponding expression c — (a + b) 
equals — ^ — m. Therefore each hypergeometric term blows up as a ^ 1. This 
divergence is made evident by employing the relation 



(12.16) 



2^1 (a, b, c; ^) = (1 - ^)'=-"-''2i='i(c -a,c-b,c; z) 



The expression for No^4{a,m) given above is transformed into 



No,i{ci,m) = 



—m — 



4r(rn+l)(l-a2)™+i/2 21^ 1 



I + 



ar(|)r(m+f) 



2T{m + 1)(1 - 02)^+1/2 I " 3 



for < 1, 



r(i)r(m + i) 



+ 



2V2ar{m + 1)(1 - a-2)'"+i/2 

r(-i)r(2m + |) 



27^1 M-"^ 1 
2i^l ( 4 



m 



2(2a)2™+3/2r(m + 1)(1 - a-2)™+i/2 
Introduce the functions 



m - 



2a[\ 



m+l 



and 



G2{a, m) 



(2a) 



2rri+l 



(-l)'"m!2-2'"(^'^+^ ) 2F1 I 4 



i 1 



m + ■ 



\ 1 



for > 1. 
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Then 

77 Gi{a,m) 



(12.17) No,4{a,m) 



4m! (1 - a2)™+V2 
for <1 and 



(12.18) iVo,4(a,TO) 



22^+5/2 ^m! (a2 - 1)^+1/2 



for > 1. The functions Gi{a,m) and G2{a,m) match at a = 1 to sufficiently 
high order to verify the continuity at a = 1. Morever, their blow up at a = —1 is a 
reflection of the fact that the convergence of the integral No,4{a, m) requires a > — 1. 

It is possible to show that both expressions (12.17) and (12.18) reduce to (12.2). 
The details will appear elsewhere. 



13. Integrals from Feynman diagrams 

The flexibility of the method of brackets is now illustrated by evaluating exam- 
ples of definite integrals appearing in the resolution of Feynman diagrams. The 
reader will find in [41], [63], [42] and [75] information about these diagrams. The 
mathematical theory behind Quantum Field Theory and in particular to the role 
of Feynman diagrams can be obtained from [34] and [22] . 

The graph G contains propagators or internal lines, L loops associated to 
independent internal momenta Q := {Qi,-- - ,Ql}-i E independent external mo- 
menta P := {Pi, • • • ,Pe} (therefore the diagram has E + \ external lines). The 
momentum Pj , Qj belong to and the space is equipped with the Minkowski 
metric. Therefore, for ^, B e M^, we have 

(13.1) := ^2 _ ^2 _ ^2 _ ^2^ 
and 

(13.2) A-B:= AqBq - A^Bi - A2B2 - A3B3. 

Finally, each propagator has a mass ruj > associated to it, collected in the vector 
m = (mi, - • • ,m7v). 

The method of dimensional regularization (see [62] for details) gives an integral 
expression in the momentum space that represents the diagram in £> = 4 — 2e 
dimensions. In Minkowski space the integral is given by 



(13.3) G = G(P,m):=|0i...| 



d^Qi 1 1 



(52 _ ^2)a, (52^ _ ^2^) 



The symbol Bj represents the momentum of the j-th propagator and it is a linear 
combination of the internal and external momenta P and Q, respectively. The 
vector a := (ai,--- ,ajv) captures the powers of the propagators and they may 
assume arbitrary values. 

In order to simplify (13.3), we use the identity 

(13.4) — = — - / x«-ie-^^ dx 
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with A = Bj — mj and convert it into 



(13.5) 




where dx = ^dxj. 



=1 



The next step in the reduction process is to integrate (13.5) with respect to 

the internal momenta Qj. This gives an expression for the integral G in terms 
of only the external momenta Pj and the masses rrij. This step can be achieved 
by introducing the Schwinger parametrizaton (see [35] and chapter 1, section 4 of 
[22] for details) and called the Schwinger variables. The final result is the 

representation 

The function F corresponds to a quadratic structure of the external momentum 
defined by 

E 

(13.7) F= ^ a,,P,-P,-. 

The function U and the coefficients Cij- are the Symanzik polynomials in the 
Sc;liwinger parameters Xj. These polynomials are given in terms of determinants 
of the so-called matrix of parameters. The polynomial Cij are symmetric, that is 
Cij = Cj^i. A systematic algorithm to write down the expression (13.6) directly 
from the Feynman diagram is presented in [35]. 

Example 13.1. Figure 2 depicts the interaction of three particles corresponding 
to the three external lines of momentum Pi, P2, P3. In this case the Schwinger 
parametrization provides the integral 




Figure 2. The triangle 
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G = 



D/2 /-oc /-oc j-oo ai-1 02-1 03-1 



// 

Jo Jo 



T{ai)T{a2)r{a3) Jo Jo Jo {xi + X2 + xa)'^/^ 

/ 2 , 2 I 2\ f C'l-Pf + 2C12P1 • P2 + C?22-P2 \ , , , 

X exp(a;iTOi +X2m2 +a;3mo)exp axiax2dx3. 

V X1 + X2 + X3 ' 

The algorithm in [35] and [36] gives the coefficients Cij as 
(13.8) Cii = xi{x2 + X3), C12 = X1X3, C22 = x^ixi + .T2). 



Conservation of momentum gives P3 = Pi + P2 and replacing the coefficients 
we obtain 

/_-|wD/2 (-oc <.oo /.oo 

G = / / / x^i-^x^^-ix""-! X 

Y{%i^{aj)Jo Jo Jo 

exp (ximf + X2mi + xamg) cxp ^;^,+4+x3 j 



(a;i + a;2 + 0:3)^/2 



dx\dx2dxz. 



To solve the Feynman diagram in Figure 2 it is required to evaluate the integral 
G as a function of the variables Pi, P2 € M^, the masses m^, the dimension D and 
the parameters Oj. 

We now describe the evaluation of the integral G in the special massless situation: 
mi = TO2 = = 0. Moreover we assume that P^ = P^ = 0. The integral to be 
evaluated is then 

1 ^\-D|2 r expf j^^^l P|) 

Gi = ^, ^'l] , / xl^-^xT-'xl^-^ : VJ dx,dx2dx,. 

r(ai)r(a2)r(a3) Jb,3 ^ (a;i + a;2 + ccs)^/^ 

The method of brackets gives 

^ ' ^ ' ^ ' m n2 ns "4 

where the brackets are given by 

Ai = (15/2 + 711+712+713 + 714), 

A2 = (ai + 7ii + 712), 

A3 = (02 + 713), 
A4 = (03 + 711 + 714). 

The solution contains no free indices: there are four sums and the linear system 
corresponding to the vanishing of the brackets eliminates all of them: 

nl ^ - ai - a2 - as, = + a2 + 03, rig = -02, 714 = +a\ + 02- 

We conclude that 

(-^\-D/2 

r<^ _ V I / p2\D/2-ai-02-03 V 

- r(ai)r(a2)r(a3)^'^^^ 

r(ai + a2 + 03 - §)r(§ ~ a2 - a3)r(Q2)r(f )r(f ~ ai - 02) 
r(Z) - ai - a2 - 03) 
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Example 13.2. The second example considers the diagram depicted in Figure 3. 
The resolution of this diagram is well-known and it appears in [12], [24] and [25]. 
The diagram contains two external lines and two internal lines (propagators) with 
the same mass m. These propagators are marked 1 and 2. 




Figure 3. The bubble 

In momentum variables, the integral representation of this diagram is given by 

d^Q 1 



^ J^n inD/^ - m^r^ {{P - Q)^ 

For the diagram considered here, we have U = X1+X2 and F = x\X2P^ ■ According 
to (13.6), the Schwinger representation is given by 

^ (-1)-^/^ 

r(ai)r(a2) 

.00 .00 ^a^-l^a,~l f X1X2 ^ ^ 

X / / ■ -y^exp (m (xi +X2)) exp ■ P \ dxidx2. 

Jo Ja (xi+X2yi^ \ X1+X2 J 

In order to generate the bracket series for G, we expand first the exponential 
function to obtain 

n3 1l) G^i^il:^ V 6..(P^rH m^Y^ f ^T^T dx,dx2 

Expanding now the term 

1 x^^x"^ 
(13-12) (^,+0:2)^/2+— = '^'''nD/2 + l-n2)^" 

with Ai = + ni — 712 + 713 + ni), and replacing in (13.11) yields 

(13.13) G = — — > 01,2,3,4 J, -Z\iA2Z\3, 

r(ai)r(a2) ^^^^^ r(f + 711-712) 



where 



Al = (-J + Til - 712 + 713 + 714), 
A2 = (ai + 711 + 713), 

A3 = (a2 + Til + 714). 
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The expression for G contains 4 indices and the vanishing of the brackets allows 
us to express all of them in terms of a single index. We will denote by Gj the 
expression for G where the index rij is free. 

The sum Gi : in this case the solution of the corresponding linear system is 

D 
2 

and the sum Gi becomes 

Gi = (-1)^^/2 



(13.14) n2 = ^ — ai — a2 — n-i, rig = — ai — ni, = —a2 — ni, 



(_^2^£)/2-ai+a2 

r(ai)r(a2) 

r(ai + 02 - D/2 + m) r(ai + ni) T{a2 + ni) [Tfi^ J 
r(oi+a2 + 2ni) nil 



oo 

X 



E 



This can be expressed as 



^ ^ ^ ' \^{ai+a2 + l), 5(01+02) 

where 

(13.16) Ai = (_i)-z./2r(ar + apD^_ 

i (oi + 02; 

The sum G2- keeping n2 as the free index gives 

nl = Y — 0-1 - 0-2 - n2, nl = a2 - Y + 1^2, nl = ai - Y + ^2, 
which leads to 

G2 = A2(Pf )^/^--+-3F2 ( " J 1 ^ „ . 

where the prefactor A2 is given by 
A2 = (-1) 



ai+a2-§, i(l + ai+a2-£>), ^(2 + oi + 02 - i:*) 

D 

2 ' - 1 ■^z 2 



4m 



p2 



.^^2 r(ai +a2- D/2)ri^ - ai)r(f - 02) 



r(ai)r(a2)r(D- 01-02) 

The cases G3 and G4 are computed by a similar procedure. The results are 
G3 = A3(Pf )-°H-mT/2-»^3F2 i(l + «i-a2), i(2 + oi-02) 

\^1 + Oi — O2, 1 — 02 + Y 

and 



4m^ 



p2 



G4 = A4(P^)--(-mT/'--3F2(', i(l-«i+«2), 1(2-01+02) 4m 

\1 — ai+02, 1 — Oi ' 



2 



p2 



where the prefactors A3 and A4 are given by 



(13.17) M = (-i)-/^S£J-^ = ,_„-./.Q£ip2Z5). 

r(a2) r(ai) 

The contributions of these four sums are now classified according to their region 
of convergence. This is determined by the parameter p := \4:m'^/P'^\. In the region 
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p > 1, only the sum Gi converges, therefore G — Gi there. In the region p < 1 the 
three remaining sums converge. Therefore, according to Rule 3.4, we have 



(13.18) G 



Gi for p > 1, 

G2 + G3 + G4 for p<l. 



We have evaluated the Fcynman diagram in Figure 3 and expressed its solution in 
terms of hypergeometric functions that correspond naturally to the two quotient of 
the two energy scales present in the diagram. 



Example 13.3. The next example shows that the method of brackets succeeds in 
the evaluation of very complicated integrals. We consider a Feynman diagram with 
four loops as shown in Figure 4. 




Figure 4. A diagram with four loops 



The methods described in [35] for the Schwinger representation (13.6) of this 
diagram, give 

U = xixsx^xr + X1X3X5XS + xixsxexr + xiX4X5Xr + X2X3X^xt+ 

XiX^XqXs, + XiX4,Xr,Xs + XiX4XQXr + X2X3X5XS + X2X3XQXr+ 

X2X4X^X'j + .T1.X3.T7.X8 + .Tia;4.T6.X8 + .T1.T5.T6X7 + T2.T3.Tga;8 + 
(13.19) X2X4X^X% + X2X4XqX'j + X^X^Xr^X-j + XiX^X-jX^ + Xi aJsXgXg + 

X2X3X'jXs, + X2X4,XeXs + X2X5XeX7 + XsXiX^Xs + X3X4X6Xr + 
T1.T5.T7T8 + T2T4T7T8 + T2T5.T6.T8 + .T3.T4.T6T8 + X3X5X6Xf + 

T2X5T7.T8 + .T3a;4X7a;8 + X3X5XQXS + X3X5XYXS 
and for the function F in (13.6): 



F = (T1T2T3T5T7 + TiT2T3a:5a;8 + a;iT2.T3.T6T7 + T1T2.T4.T5T7+ 
0:1X2X3X6X8 + X1X2X4X5X8 + X1X2X4X6X7 + X1X3X4X5X7+ 
X1X2X3X7X8 + X1X2X4X6X8 + X1X2X5X6X7 + X1X3X4X5X8+ 
X1X3X4T6X7 + X1X2X4X7X8 + T1X2.T5T6X8 + X1T3T4X6X8+ 
X1X3X5X6X7 + X1X2X5X7X8 + X1X3X4X7X8 + X1X3X5X6X8+ 
X1X3X5X7X8) P^. 
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The large number of terms appearing in the expressions for U and F (34 and 21 
respcetively) makes it almost impossible to apply the method of brackets without 
an apriori factorization of these polynomials. This factorizations minimizes the 
number of sums and maximizes the number of brackets. In this example, the 
optimal factorization is given by 

F = ari/7 

(13.21) 

U = xife + fr, 
where the functions fi are given by: 

fe = xzji + (2:4/4 + /s), 

/5 = X3(x4/4 + h), 

(13.22) /4 = X5/2 + (X6/2 + /i), 

h = xz{x&h + /i), 

/2 = {X-J + Xs), 
/l = X7Xi. 

To analyze the diagram considered here, we start with the parametric represen- 
tation 



( v\-D/2 7 exp f , f Pi 

(13.23) G = ^-^1— / V ^ih + h 

n r(a,) (^1-^6 + Ii) ' 

and expand the exponential function first. A systematic expansion associated to 
the polynomials fi leads to the order 



(13.24) u h-^ U-^ h-^ /2, 

that yields the bracket series 



(13.25) G-i^J '/'ni>->".aT;7^^f^wn^' 

where we have defined the factor 



= . 

r(-ni - n3)r(-n2 - n4)r(-n5 - n7)T{-n% - n8)T{-ng - nii)r(-nio - ni2) ' 

and the corresponding brackets by 
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(13.26) 



Ai 


= + ni + n2 + na) , 


Ag = { 


[a2 4 


-jn) , 




A2 


= (-ni - na + 714 + "5) , 


Aio = 


(as - 


f ns -f 




A3 


— (— «2 — n4 + ne + n-i) , 


All = 


(a4 - 


-f ns) , 




A4 


= (-^^5 -nr + ns + ng) , 


A12 - 


(as - 


f ng 4- 


- nio) , 


As 


= {-ne -ns + nio + nn) , 


Ai3 = 


(ae - 


■f ni2) 


7 


As 


= (-ng - nil + ni2 + nu) , 


Ai4 = 


{ar - 


■f ni3 - 


-f ni4) , 


A7 


= (-nio - ni2 + ni4 + nis) , 


Ai5 = 


{as - 


f ni3 - 




As 


= (ai + ni + n2) . 











There is a unique way to evaluate the series: the numbers of indices is the same as 
the number of brackets. Solving the corresponding linear system leads to 



,p.>„; .nr(-»5) 



where the values n* are given by 



n r(a,) 



Til = 2D — ai — a2 — as — Ui — — uq — a-j — as, 
nl = a2 + aa + a4 + a5 + ae + ar + as - 2D, 
n*3 = ai- f , 
nl = -02, 

ng = ^ - as - 04 - 05 - ae - a? - as, 
ng = a4 + ag + ae + 07 + as - 

= aa- -j, 
nt = —04. 



as - 06 - ar - as. 



n*o=D 



nio = ae + a7 + as - £>, 

<i = ag - f , 



n: 



-ae, 



nl3 = 



(■14 



07 - 
n D 



as. 



Example 13.4. The last example discussed in this paper gives the value of a Feyn- 

man diagram as a hypcrgcomctric function of two variables. The diagram shown 
in Figure 3 contains two external lines and with internal lines (propagators) with 
distinct masses. The same diagram with equal masses was described in Example 
13.2. The integral representation of this diagram in the momentum space is given 

by 



(13.28) = /iJ 



1 



/2 (Q2 - m2)«i {{P - Q)2 - m2^"^ ' 



2J 



On the other hand, the parametric representation of Schwingcr is 

nr(a,)/ (^1 + ^2)- 



In order to find the bracket series associated to this integral, we first expand the 
exponentials 
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r ni+n3 n2+n3 

"1 2^ ^m.n^ns (-TOl) (-TO2) (-P ) / ; -n^^ 

n r(aj) "i."2,"3 (a;i + a;2) ^ 



-'J 

and then the denominator 



1 ^"4„n6 

We obtain 

n r(aj) "1'-'"^ -LU+'^sj 
Using the rules for transforming integrals into brackets yields 

(13.30) <^-^ 2^ 0ni,...n5 p/D , ^ X A1A2A3, 

nr(a,)-.-."a i(2+«3) 

with brackets defined by 

Ai = + ns + n4 + ns) , 

(13.31) A2 = (ai + ni + ns + 714) , 

A3 = (a2 + + ns + "5) • 

We have to choose two free values from ni, • • • ,715. The result is a hypergeometric 
function of multiplicity two. There are 10 such choices and the brackets in (13.31) 
produce the linear system 

= § + ns + n4 + n5 

(13.32) = ai + ni + ng + 714 

= 02 + 712 + 7^3 + n5. 

We denote by Gij the solution of this system with free indices ti, and Uj. The 
series appearing in the solution to the Feynman diagram in Figure 4 is expressed 
in terms of the Appell function F4, defined by 



(13.33) 









(; 








s 





„ {j)m{S)n m\ n! ' 

m.n— 



E 



Following a procedure similar to the one described in the previous example, we 
obtain the explicit values of the integral G is given in terms of the functions Gij . 
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G 



Gi,2 + Gi,4 + G2,5 when mf, ml < 
Gi,3 + Gs,5 when , < ml 

^2,3 + G3,4 when ml,P'^<mf. 

These in turn are expressed in terms of the Appell function: 



Gi,: 



(-1) 



^ r(ai+a2-f)r(f -ai)r(f 

r(ai)r(a2)r(i) - oi - 02) 



1 + ai + a2 — D 



l + oi - 



ai + 02 — 
1 +02 - 



2 2 

p2 ' p2 



r(a2-# 



:f2,5 




G2 • 



D r(ai - f ) , 2\f -<il / 2\ 



r(ai) 



2 



2 



02 



1 - 01 + ■ 



mf 



.£r(ai+a2-f)r(^ 



■ 02) 



r(ai)r(f] 



-mi) 



ai + 02 — 



a2 



2 



ma 



-#r(a2-f) 



r(a2) 



1 + 02- 



/ 2\ -s-""^ c 2\ - 



D 
2 



ai 



? 1 - a2 + ? 



-)2 2 

ma 
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14. Conclusions and futtjre work 

The method of brackets provides a very effective procedure to evaluate definite 
integrals over the interval [0,oo). The method is based on a heuristic list of rules 
on the bracket series associated to such integrals. In particular we have provided a 
variety of examples that illustrate the power of this method. A rigorous validation 
of these rules as well as a systematic study of integrals from Feynman diagrams is 
in progress. 

Acknowledgments. The authors wish to thank R. Crandall for discussions on an 
earlier version of the paper. 
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